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We compute the hadronic matrix elements of the four-quark operators needed for the study 
of K° — K° mixing beyond the Standard Model (SM). We use nf = 2 + 1 flavours of domain-wall 
fermion (DWF) which exhibit good chiral-flavour symmetry. The renormalization is performed non- 
perturbatively through the Rl-MOM scheme and our results are converted perturbatively to MS. 
The computation is performed on a single lattice spacing a ~ 0.086 fm with a lightest unitary pion 
mass of 290 MeV. The various systematic errors, including the discretisation effects, are estimated 
and discussed. Our results confirm a previous quenched study, where large ratios of non-SM to SM 
matrix elements were obtained. 



Introduction. — Recent progress achieved by the lat- 
tice community is greatly improving our theoretical un- 
derstanding of CP violation in kaon decays. The exper- 
imentally well- measured parameters Ek and e' Kl which 
quantify indirect and direct CP violation in K — > tttt 
can be confronted with a theoretical computation of the 
K — > (7T7r)i = o.2 amplitudes and of neutral kaon mixing, 
providing the non-perturbative effects are correctly ac- 
counted for. The first direct and realistic computation 
of the K — > (7T7r)i = 2 amplitude has been recently per- 
formed [HE]- A complete computation of the AI =1/2 
amplitude is still missing, but important work is being 
performed in that direction [3]. The situation is much 
more favourable for neutral kaon oscillations: the bag 
parameter Bk which describes the long-distance contri- 
butions to neutral kaon mixing in the SM, is now com- 
puted with a precision of a few percent (see e.g. [1H^) for 
recent unquenched determinations). By combining 
with the experimental value of sr, one obtains impor- 
tant constraints on the free parameters associated with 
quark flavour mixing (see [7] for a recent pedagogical re- 
view). In principle, the same techniques can be applied 
for beyond the Standard Model (BSM) theories (see for 
example [5HPT]). but not much is known concerning the 
long-distance contributions of the non-standard opera- 
tors beyond the quenched approximation. 

In this letter, we present the first realistic compu- 
tation of the matrix elements of neutral kaon mixing 
beyond the standard model. Previous studies were 
either preliminary [181 119] or suffered from the quenched 
approximation [201 [21] • Since a noticeable disagreement 
was observed between [20] and [21] it is important to 
repeat this computation in a more realistic framework 



Formalism. — In the SM, neutral kaon mixing is dom- 
inated by box diagrams as in figure [T] By perform- 
ing an operator product expansion, one factorises the 
long-distance effects in the weak matrix element (WME) 
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FIG. 1: Diagram contributing to A' — A' mixing in the SM. 



(K°\O^ s ~ 2 \K ) of the four-quark operator Of s_2 given 
by (a and /? are colour indices) 

Of 8=2 = (*a7„(l ~ lM (Spj^l ~ 75 )cW • (1) 

Only one four-quark operator appears in the SM because 
neutral kaon mixing occurs under W-boson exchange, im- 
plying a "vector-axial" Dirac structure. Since this four- 
quark operator is invariant under Fierz re-arrangement, 
the two different colour structures (mixed and unmixed) 
are equivalent. 

We also seek to understand whether new physics be- 
yond the standard model could play a detectable role 
in kaon CP violation. Assuming that this new physics is 
perturbative, we might also describe this by contributions 
to the effective Hamiltonian, and other four quark oper- 
ators are induced by such extensions of the SM. It is con- 
ventional to introduce the so-called SUSY basis 0^?= 2 5 : 
in addition to the SM operator Of 3=2 , we define [22] jll] 
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These four-quark operators appear in the generic As = 2 
Hamiltonian 



#A S =2 = ^ C . (M)0 A s =2 (/i) 



(3) 
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where /i is a renormalization scale. The Wilson coeffi- 
cients Ci, which encode the short-distance effects, depend 
on the new physics model under consideration. The long- 
distance effects are factorised into the matrix elements of 
the four-quark operators 0^ s=2 . Lattice QCD offers a 
unique opportunity to quantify these effects in a model- 
independent way. 

In phenomenological applications, by combining our 
results for the BSM matrix elements with experimen- 
tal observables (typically the mass difference AMjf = 
rriK L — ttik s and ck), one obtains important constraints 
on the model under consideration. The Wilson coeffi- 
cients and the bare matrix elements have to be converted 
into a common scheme, at a given scale \x (typically MS 
at a scale of 2 or 3 GeV). In the SM case (i = 1) it is 
conventional to introduce the kaon bag parameter Bk, 
which measures the deviation of the SM matrix clement 
from its vacuum saturation approximation (VSA) 



with an inverse lattice spacing a -1 = 2.28(3) GeV, corre- 
sponding to a ~ 0.086 fm. The quarks are described by 
the domain- wall action 25 27J, both in the valence and in 
the sea sectors. We have three different values of the light 
sea quark mass am^ ht — 0.004,0.006,0.008 correspond- 
ing to unitary pion masses of approximately 290, 340 and 
390 MeV respectively. The simulated strange sea quark 
mass is an^tfangc = 0.03, which is close to its physical 
value 0.0273(7) [21] ■ For the main results of this work, we 
consider only unitary light quarks, am^ ; 
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tially quenched (am™ 1 ™ = 0.025) data. 

To extract the bare matrix elements we follow [H [21] , 
where the procedure for the evaluation of the two-point 
functions and of the three point function of the SM op- 
erator has been explained in detail. In particular, we 
have used Coulomb gauge fixed wall sources to obtain 
very good statistical precision. From the measurement of 
Bk, the generalisation to the BSM operators is straight- 
forward. We define the three point function 



Where the normalisation for the decay constant is such 
that fx- = 156.1MeV. Several normalisations for the 
BSM operators have been proposed in the literature (see 
for example [10]), in this work we follow [ST] and define 
the ratios 
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where P is a pseudo-scalar particle of mass m p and decay 



constant f P . The term [l Iatt is obtained from our lattice n S htest unitar y ^ n am ™i* 
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simulations for different values of mp. 

As advocated in [21] . there are various reasons to 
choose this normalisation rather than, for example, the 
standard bag parameters. In particular, we expect some 
systematic errors to cancel in the ratio of the bare WMEs, 
there is no need to introduce the light quark masses and 
the partial conservation of the axial current. Further- 
more, the mass factors in eq.(|5| compensate the leading- 
order chiral behaviour of the WMEs, making all of the 
chiral extrapolations smoother. Finally, at the physical 
point P = K° , the Rf SM 's give directly the ratios of the 
non-SM to SM contributions. 

The Rf s will be the main result of this work, but for 
completeness we will also give the BSM bag parametri- 
sation, defined as (where N2,,..,s = §,— §,— 2, — f) [23] . 



From the asymptotic Euclidean time behaviour of 
the ratios of three point-functions Cj/ci (we fit these 
ratios to a constant in the interval t/a = [12,52] 
[45] ) we obtain the bare matrix elements of the 
four-quark operators normalised by the SM one: 
[(P\0^ s = 2 \P)/(P\0^ s=2 \P)] harc . In figure § we 
show the corresponding plateaux obtained for our 
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FIG. 2: Ratios of the bare three point functions from which 



Computation details. — This computation is performed 
on the finer of the two ensembles described in detail 
in [U[2D. ^ e use ^ 3 x 64 x f6 (the last number corre- 
sponds to the length of the fifth dimension of the domain- 
wall action in lattice units) Iwasaki gauge configurations 



we extract Ri 
unitary kaon. 



Results are shown for our lightest simulated 



Renormalization. — The four-quark operators given in 
eq. (f|2) mix under renormalization. In a scheme which 
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preserves chiral symmetry, the mixing pattern is given 
by the transformation properties of these operators un- 
der chiral rotations SU(3) L x SU(3) R . The SM opera- 
tor 0^ s=2 belongs to a (27, 1) irreducible representation 
of SU(3)l x SU(3)r and renormalizes multiplicatively. 
The BSM operators fall into two categories: O^" 2 and 
2 transform like (6, 6) and mix together. Likewise 
2 and 0^ s=2 belong to (8, 8) and also mix. If chi- 
ral symmetry is realised, the five-by-five renormalization 
matrix Z^ is block diagonal: a single factor for the (27, 1) 
operator and two 2x2 matrices for the BSM operators. 
Because we work with the domain-wall fermions formu- 
lation, in which the explicit breaking of chiral symmetry 
can be made arbitrary small (and in practice numerically 
irrelevant) , we expect to find this continuum pattern, up 
to small discretisation effects. 

We perform the renormalization of the four-quark 
operators 0^ s=2 non-perturbatively in the RI-MOM 
scheme [2HJ- We compute the forward, amputated, ver- 
tex functions of the relevant operators between external 
quark states in the lattice Landau gauge for a given set 
of momenta. As a renormalization condition, we im- 
pose these to be equal to their tree-level values once 
projected onto their colour-Dirac structures and extrap- 
olated to the chiral limit. By using both momentum 
sources [29] and partially twisted boundary conditions, 
we obtain smooth functions of the external momentum 
with very good statistical accuracy 05] ■ Although we 
have performed this computation in various intermediate 
schemes (including the non-exceptional (7^, 7^ ^scheme- 
introduced in [4]) we quote here the results obtained via 
the RI-MOM scheme because only in this case are the 
conversion factors to MS (computed in continuum per- 
turbation theory) available for the whole set of operators. 
In such a scheme, the presence of exceptional channels en- 
hances the Goldstone pole contaminations p~8| . I30 j . which 
have then been subtracted explicitly [31]. We choose 
to impose the renormalization conditions at \i = 3 GeV. 
At this scale, perturbation theory (PT) converges rather 
quickly (a s ~ 0.25), the chiral symmetry breaking ef- 
fects are small, and we still have good control on the 
discretisation effects [4j [32] . Matching at this scale is 
important for our error budget since the perturbative 
conversion factors between the lattice schemes and MS 
of the four-quark operators are only known at one-loop 
(as discussed later, matching to PT is actually one of 
our dominant sources of systematic error). We observe 
that the effects of chiral symmetry breaking are not com- 
pletely negligible, even at 3 GeV [33J . Therefore we must 
assess a systematic error to the mixing of operators of dif- 
ferent chirality (see next section). We have checked that 
in a non-exceptional scheme this small chirally forbidden 
mixing is strongly reduced and becomes numerically ir- 
relevant at 3 GeV [TH1 [33J . Thus we conclude that this ef- 
fect is a physical manifestation of the infrared behaviour 
of the exceptional intermediate scheme. The results for 
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FIG. 3: Renormalization factors of the four-quark operators. 
At each value of the simulated momentum p, we run to the 
scale of 3GeV and convert to MS. The remaining scale depen- 
dence can be imputed to the truncation of the perturbative 
expansion. We show only the ^-factors allowed by chiral sym- 
metry. 



the chirally allowed renormalization factors Z^ s (3GeV) 
are shown in figure [3] They relate the bare four-quark 
operators to the renormalized ones through the usual re- 
lation (Z q is the renormalization factor of the quark wave 
function and cancels in the ratios) 
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Physical results and error estimation. — Once the bare 
ratios have been renormalized, we extrapolate them to 
the physical kaon mass. We have explored different 
strategies, such as a simple polynomial form or next-to- 
leading order chiral perturbation theory (ChPT) predic- 
tions [47 . The value of the simulated unitary strange 
mass differs somewhat from the physical one, therefore 
we perform an interpolation using a partially quenched 
strange quark. We find that the i? BSM 's exhibit a very 
mild quark mass dependence (see figure [4] [48]), therefore 
we take the results obtained with the polynomial Ansatz 
as our central values [3FJ. 

Our final results are the Rf SM quoted in MS at 3 GeV 
given in table [I] For completeness, we also convert these 
to the BSM bag parameters, using eq. Q 50J . We 
also note that, using the same framework, the SM con- 
tribution is found to be B\ = Bk = 0.517(4) sta t in 
the MS scheme at 3 GeV, whereas a continuum value 
of 0.529(5) s tat(19) S yst was quoted in 0]. The difference 
comes from the fact that a different intermediate scheme 
was used in [4] (such a difference is accounted for in our 
estimation of the systematic errors). From the same ref- 
erence, the discretisation effects for Bk on this lattice are 
seen to be of the order of 1.5%. Since we have only one 
lattice spacing for the BSM ratios, we make the assump- 
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FIG. 4: Renormalized BSM ratios -Rf SM in function of the 
bare valence quark mass. We show the three unitary light 
quarks for both the unitary strange (upward blue triangles) 
and the partially quenched strange (downward blue triangles), 
together with their extrapolations in the light sector (blue and 
red circles) and their interpolation to the physical kaon mass 
(black squares). 



tion that the discretisation errors are of the same size 
as those affecting Bk, and estimate a 1.5% error to all 
the different operators. This might appear like a crude 
estimate, but this effect is expected to be sub-dominant 
compared to other sources of systematic errors. The next 
systematic error (called "cxtr.") represents the spread of 
the results obtained from different extrapolation strate- 
gies to the physical point. The systematic error associ- 
ated with the non-perturbative renormalization (NPR) 
has been estimated from the breaking of chiral symme- 
try. The mixing of the (6, 6) with the (8, 8) operators is 
forbidden by chiral symmetry, but likely to be enhanced 
by the exchange of light pseudo-scalar particles. As the 
matrix element of 0^ s=2 is numerically large, this non- 
physical mixing has effect on 0^ s=2 and 0^ s=2 of the or- 
der of 8 — 9%. As discussed in the previous section, this 
unwanted infrared effect is absent if a non-exceptional 
scheme is used. The last error we quote ("PT") arises 
from the matching between the intermediate RI-MOM 
scheme and MS, which is performed at one- loop order in 
perturbation theory [34, 35 in the three-flavour theory. 
The associated error is obtained by taking half the differ- 
ence between the leading order and the next to leading 
order result [5T]. We note that this error is one of the 
dominant ones in our budget, and we expect this error to 
be reduced by an important factor if a non-exceptional 
scheme were used, since the latter are known to converge 
faster in perturbation theory. We neglect the finite vol- 
ume effects which have been found to be small in as 
one can expect from the value of m^L ~ 4 for our lightest 



i 


i?f SM Bi 


stat. 


discr. 


extr. 


NPR 


PT 


total 


2 


-15.3(1.7) 0.43 (5) 


1.3 


1.5 


4.0 


9.4 


4.7 


11.3 


3 


5.4(0.6) 0.75 (9) 


2.0 


1.5 


3.9 


7.8 


7.6 


12.0 


4 


29.3(2.9) 0.69 (7) 


1.3 


1.5 


4.1 


3.0 


8.2 


9.8 


5 


6.6(0.9) 0.47 (6) 


2.1 


1.5 


3.8 


3.2 


12.6 


13.8 



TABLE I: Final results of this work: the first two columns 
show the ratios Rf SM and the corresponding bag parameters 
Bi in MS at 3 GeV, together with their total error, combining 
systematics and statistics. In the remaining columns, we give 
our error budget for the R BSM , detailing the contributions in 
percentage of the different sources of systematics (see text for 
more details). 

pion mass ~ 290 MeV. 

Conclusions. — We have computed the electroweak 
matrix elements of the As = 2 four-quark operators 
which contribute to neutral kaon mixing beyond the SM. 
Our work improves on the the previous studies by using 
both dynamical and chiral fermions. We confirm the ef- 
fect seen in a previous quenched computation |21j . where 
a large enhancement of the non-standard matrix elements 
were observed. The errors quoted in this work are of 
the order of 10%. We note that the main limitation of 
this study comes from the lack of matching factors be- 
tween non-exceptional renormalization schemes (such as 
SMOM) and MS. Once these factors are available, we 
expect to reach a precision better than 5%. We also plan 
to utilise another lattice spacing in order to have a better 
handle on the discretisation effects. 



Appendix 

We have computed the non-perturbative scale evolu- 
tion of the i? BSM 's between 3 and 2 GeV, and then con- 
verted the results to MS using one-loop perturbation the- 
ory [MUSS]: 
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Our conventions are such that 

i? BSM (2GeV) = U m (2GeV,3GeV)R BSM 



(3 GeV) . 

(10) 
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